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ABSTRACT 

Using high-resolution N-body simulations, we examine whether a major dry 
merger mitigates the difference in the radial density distributions between red 
and blue globular clusters (GCs). To this end, we study the relation between the 
density slope of the GCs in merger progenitors and that in a merger remnant, 
when the density distribution is described by uqq oc r~ a . We also study how 
our results depend on the merger orbit and the size of the core radius of the 
initial GC density distribution. We find that a major dry merger makes the GC 
profile flatter, and the steeper initial GC profile leads to more significant flatten- 
ing, especially if the initial slope is steeper than a ~ 3.5. Our result suggests 
that if there is a major dry merger of elliptical galaxies whose red GCs have a 
steeper radial profile than the blue GCs, as currently observed, and their slopes 
are steeper than a ~ 3.5, the difference in the slopes between two populations 
becomes smaller after dry mergers. Therefore, the observed slopes of red and 
blue GCs can be a diagnostic of the importance of dry merger. The current 
observational data show that the red and blue GCs have more comparable and 
shallower slopes in some luminous galaxies, which may indicate that they have 
experienced dry mergers. 
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Introduction 



Globular clusters (GCs) in elliptical galaxies have been intensively studied in consider- 
atio n of explaining the form ation of both their host elliptical galaxies and GCs themselves 



see 



Brodie fc Straderl 12006k for a review ). GCs are attractive as tracers of the star for- 



mation history of their host galaxies (e.g. lYi et al.l 12004 ; IStrader et al.ll2006l ). because some 
properties of GC systems are correlated with the p roperties of their host galaxies (e.g., 



Brodie fc Huchral [l99l| ; iDjorgovski fc Santiago! Il992h . It is thought that t he formation of 



Ashman & Zepi 


1992) 


(Forbes et al. 


1997 


)• I 



star clusters like current old GCs in local galaxy mergers (ISchweizerl 120061 ) . 



An important aspect of GC systems in elliptical galaxies is a color bimodality (e.g 



Zepf fc Ashman 



Peng et al 



1993 



Geisler et al.lll996l : iGebhardt fc Kissler-Patidll999l : lLarsen et al.ll2001 



2006 1). It is also found that red GCs are more centrally c oncentrated than blue 



GCs (e.g. iForte et al.ll2005l ; iBassino et al.ll2006al : iTamura et a 



20061) . which could put ad- 
ditional constraints on their formation scenario (jBekki et al.ll2002l ). The radial profile of 
each GC subpopulation is well described by a power-law distribution, especially at the outer 
radii, and the red GCs have a steeper slope than the blue GCs. Although the origin of 



Yoon et al. 


2006: 


Strader et al. 


2007 





20071 ). Classically, three scenarios have been proposed to explain the color bimodality; ma- 



jor gas-rich mergers, in situ formation of multiple dissipat i onal c ollapses, and dissipationless 



accretion. An explanation suggested by lAshman &: Zepj (119921 ) is that red GCs are metal- 



rich clusters which might be formed by gas-ric h disk-disk mergers. Therefore, the red GCs 



might be younger than blue GCs. Meanwhile, iForbes et al.l (119971 ) explain that blue GCs 



might be formed in the first stage of dissipational collapse and red GCs might be formed 



after the truncation of the blue GC formation. Another explanation given by ICote et al. 



( 119981 ) includes accretion of b l ue GC s from small galaxies into the already formed red GCs. 
More recently, iBeasley et al.l (120021 ) demonstrate that the bimodality can be explained in 
elliptical galaxy formation based on a hierarchical clustering scenario. 



On the other hand, recent research suggests that in the late stage of evolution, early-type 
galaxies might have experienced dry merging where merger progenitors do not have much 
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gas, nor accompany star formation. The number density evolution of red galax ies has been 



discu ssed in observations of COMBO-17 (jBell et al.ll2004l ) and DEEP2 surveys (IFaber et al. 



20051 ). and such studies suggest that th e density change can be understood by the dominance 
of dry merging after z < 1 (but see also Yamada et all200~5 ; Cimatti et all2006 ; Bundy et al. 



20071 ; IScarlata et al.ll2007l ). The evolution of galaxy clustering a 



evolution of massive red galaxies from dry merging (IWhite et a 
yations show that dry merging does occur Jvan Dokkumll2005l : 



so im plies late effects on the 



2007). Moreover, the obser 



Tran et al. 2005; Rines et al. 



20071), while the observe d features of galaxies are well explained in cosmological simulations 



(e.g. iKawata et al.l 120061 ). Recent theoretical studies of dry merger simulations of ellipticals 
show that merger re mnants maintain their properties on the fundamental plane and other 
scaling relations 



e.g. 



Nipoti et al.ll2003l : iBovlan-Kolchin et al.ll2005l : iRobertson et al.l 12006 



Ciotti et al.l 120071 ). T he dry merging o f binary ellipticals also can explain the formation of 
boxy- type ellipticals (INaab et al.l 120061 ). 



Bekki fc Forbes! (120061 ) demonstrate that the observed correlation between a spatial 



distribution of GCs and the total luminosity of ellipticals can be explained by sequential 
dissipationless major mergers, because the radial density profile of GCs progressively flatten 
after each major dry merger. This study raises the important question of whether or not 
the slopes of the density profiles of red and blue GCs persist after major dry merging. For 
example, we now consider that the density profiles of red and blue GCs in progenitor elliptical 
galaxies are described by n GC red oc r~ a P' rod and n GC b i ue oc r _a p> bluc , and these profiles in a 
major dry merger remnant become necked oc r~ ai ' rcd and riGC,biue oc r _a ^ blue . The current 
observations suggest that a r)re( j > a r ,biue- However, if a dry merger flattens the red GC 
density profile more than the blue GC density profile, a Pjre d — a rjre d > a Pi bi uc — ^r.biue, in 
the remnant galaxy the difference between a r)re d and a rj biue becomes smaller. In this case, 
the observed difference in the slopes of the density profiles of red and blue GCs in nearby 
ellipticals can be a valuable diagnostic for the importance of the dry merger in the evolution 
of ellipticals. To clarify this issue, the question becomes how the flattening of GC profiles 
during dry merging depends on the initial distributions of the GCs. 

We use numerical simulations of major dry mergers to study the dependence of GC 
distributions in merger remnants on the initial distributions in merger progenitors. Then, 
we can compare a Pire d — a riI . e d and a Pi biue — o^biue, for different sets of a Pjre d and a Pjb i ue . 
Since major dry mergers between two equal-mass merger progenitors must leave the most 
significant effects on merger remnants, compared with minor mergers, we study only equal- 
mass mergers in this paper. In §2, we explain details of our dry merger simulations and 
initial GC distributions. The changes in GC spatial distributions due to major mergers are 
shown in §3. We discuss the implication of our results in §4 that is followed by conclusion. 
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2. Simulation and parameters 



Our galaxy merger simulations are carried out with a Tree N-body code, GCD+ ( IKawata fc Gibson 



20031 ). We model a merger progenitor elliptical galaxy which consists of a stellar bulge and 
a dark matter (DM) halo. The stellar component is assumed to follow the iHernquistl (119901 ) 
density profile, 

n M * a 1 m 

P * W = 2^r (7T^ W 

10 11 M(T). Following iBoylan-Kolchin et al.l (120061 ). we use 



We set the total stellar mass M* 



an empirical relation from IShen et al.l (120031 ) , 



R P = 4.16 



0.56 



10 n M Q 

For the Hernquist profile, R e can be described with R e 



kpc. 



(2) 



1.8153 a. 



The initial DM halo is assumed to follow the NFW density profile (jNavarro et al.lll9971 ); 

Q 5 C 



P 



8ttG 



fl(z) cx(l + cx) 2 ' 



(3) 



where c is a concentration parameter, and x = r/r 2 oo is a radius normalized by r 2 oo which is 
the radius of a sphere whose mean interior density is 200 p crit . The characteristic density S c 
is described as 



200 



3 ln(l + c) - c/(l + c) 
The radius of r 2 oo is linked with M 20 o, the mass within the radius, with 



r 200 = 1.63 x 10" 



M 



200 



1/3 



-1/3 



;i + z )- 1 /i- 1 kpc. 



(4) 



(5) 



We assume c = 10, Zo — 0; Hq — 70 km s _1 Mpc -1 . The final DM density profile takes int o 
account a diabatic contractio n , follo wing the formula suggested by iBlumenthal et al.l (119861 ). 
Following ISpringel fc White! (119991 ). we apply an exponential cut-off for the dark matter 
density profile at r > r 2 oo- We define the total mass of the DM component as the DM mass 
within 2r 2 oo, where the density is low enough. The total mass of the DM component is set 
to be M DM ,tot = M NF w(2r 2 oo) — M*, where M NF w(2^ 2 oo) is calculated from the exponentially 
truncated NFW profile with M 2 qo = 10 12 M Q . The DM and stellar components are composed 



of 10 5 and 10 4 particles, respectively, so that the particles for the different components have 
roughly same mass. As a result, the particle mass and softening length are ~ 10 7 M and 
0.58 kpc, respectively. We also run one simulation with half number of particles in order to 
verify convergence of our results. 
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The velocity dispersions for both DM and stellar components are assumed to be isotropic, 
o> = oq = <70, and follow the Jeans equation. Then, the velocity dispersion can be written 

as 

f°° d<? 
p(r)o*(r)=G p(s)M(<s)- (6) 



as 



shown in iMamon fc Lokasl J2005h . In our simulations, the two merger progenitors have 



exactly the same distributions of DM and star particles for simplicity. 

Three different merger orbits are studied in our simulations as given in Table 1 . From 



a fixed initial separation of 400 kpc, we have chosen three parabolic orbits (lBarneslll992l ); 
(A) head-on collision; (B) a close encounter with pericentric distance of the orbit, r per i = 5 
kpc; and (C) a wider encounter with r peri = 50 kpc. All simulations are stopped after 8 Gyrs 
even though the merger remnants are fully relaxed within 300 kpc within 6 Gyrs. 

For simplicity, we assume that the GCs are spherically distributed and that the number 
density follows the following radial profile: 

™Gc(r) oc (r2 + ^ 2)Qp/2 , (7) 
where r c is a core radius. We in clude the core, bec ause the observed GC density profiles 



suggest the existence of the core (IForbes et al.lll996l ). We also examine how the initial core 



size affects the slope at the outer radii in a dry merger remnant. Here, instead of putting 
additional particles to trace the positions of GCs during the merger, we randomly select 
particles from DM particles, such that the selected GC particles follow the number density 
distribution of equation ((Tj) as well as velocity dispersion distribution of equation (jSj). This 
sampling of GC particles has the advantage that we can simulate various different initial 
distributions without running a simulation several times, which enables us to study the 
detailed dependence of the flattening on the initial GC profiles. 

After the merger, we analyze the number density profile of the selected GC particles. 
We are interested in how final slope of the GC density profile depends on the initial slope. 
We fit the profile of the number density at the radii between 10 and 100 kpc by hqq oc r~ a \ 
This radial region is chosen arbitrarily. However, the important point is that we compare 
the slope at a fixed radial region throughout the paper. We consider initial density slopes 
between a p = 2 and 5, and the two different core radii of r c = 2.5 and 5 kpc to see the 
effect of the initial core size on the final slope. For each initial GC profile, we sample 10,000 
GC particles for each merger progenitor. Figure [1] demonstrates that our sampling well 
reproduces the input number density profile for GCs even at the outer radii. It also confirms 
that the distribution of velocity dispersions for the GCs is consistent with the equilibrium 
values from equation ( [6]) . 
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We also test stability of the initial distribution of the selected GC particles by running 
a simulation of a single isolated galaxy with the sampled GC particles. After 8 Gyrs, the 
isolated single galaxy maintains the initial distribution of the GCs, as shown in Figure [2J 



Results 



We present the projected density distribution of the stellar component of a merger 
remnant in Figure [3j We can see that there is a negligible effect from the difference in orbits 
on the distribution of stars in the merger remnant. We also fit the final stellar density profiles 
between r = 1 and 40 kpc with de Vaucouleurs profiles. In Figure [31 de Vaucouleurs profile 
describes well the projected distribution of the stellar component in all merger remnants 
from the three different orbits. The measured effective radius is close to ~ 6.2 kpc for all 



three simulations, which is similar to the results in iBoylan-Kolchin et al.l (120051 ). 

In Figure H] we show the final GC number density profile for the different models with 
a p = 2, 4 and 5 for simulation B. We first focus on the results for simulation B, and discuss 
the orbital dependence later. For models of a p = 4 and 5, the GC profile in the merger 
remnants is clearly flatter than the initial one. This result is c onsistent with wha t have been 



found in s i milar dissipationless galaxy merger simulations ( White! Il978l . Il980l ; IVillumsen 



20061 ). On the other hand, for 



19821 . Il983l : iDuncan et al.lll983l : IWhitelll983l : iBekki fc Forbes 
the a p = 2 model there is not much change in the density slope of the GC distribution in a 
merger remnant. In addition, comparing models of a p = 4 and 5, the difference between the 
initial and the final profile is more significant for the a p = 5 model. Hence, flattening of the 
GC radial profile due to a major dry merger clearly depends on the GC density profile in a 
progenitor galaxy. 

Figure [5] summarizes our main results. We find that a higher a p leads to a larger a p — a T . 
In other words, the steeper initial GC profiles result in stronger flattening. The boundary 
between a significant and a weak flattening is around a p ~ 3.5. In the next sections, we 
discuss the implication of our results for the evolution in the relative distributions of red 
and blue GCs during a major dry merger. First, we briefly examine the importance of the 
merger orbits and the initial core radii. 



3.1. Merger Orbits 



We find that the final GC density profile weakly depends on a merger orbit as shown in 
Figure El The found maximum difference of a r is ~ 0.2 for a p = 5. As also highlighted in Fig- 



- 7- 



ure[7]of the three orbits, the simulation C causes the least flattening in the GC distribution. 
The highest angular momentum collision is not efficient in flattening GC distributions. 

We also find that the density profile in the outer regions of merger remnants are more 
sensitive to the merger orbits than that in the inner region. In Figure [7J all three merger 
orbits do not produce any differences in r < 10 kpc for an initial a p = 5 model. However, 
beyond r ~ 10 kpc, the small difference of flattening begins to appear. 

As a convergence test, we also run simulation A with a half number of particles, i.e, 
5 x 10 4 DM, 5 x 10 3 stellar, and 5 x 10 3 GC particles. Figure [2] demonstrates that the 
low resolution simulation also provides a stable initial condition. As shown in Figure the 
low-resolution simulation leads to the consistent results to the high-resolution simulation. 
Thus, our results are not affected by the number of particles adopted. 

3.2. Dependence on the size of core radius 

In Figure [91 we compare the GC slopes after major dry mergers in the cases where the 
merger progenitors have the core radii of r c = 2.5 and 5 kpc. Larger initial core radii result 
in stronger flattening. Therefore, it is clear that the flattening due to a dry merger depends 
on the initial core radius. 

The size of the initial core radius is a more important factor for the final density slope 
than the merger orbit, as found in the comparison between Figure 161 and [TUl In particular, 
the size of the initial core radius affects the steepening of the GC density profile differently, 
depending on a p . The difference in a r between the different initial core radii is more signif- 
icant for a larger a p model than for a smaller a p . Hence it is another interesting issue to 
see in more detail how the initial core radius affects the flattening. However, in our method 
described in §2, it is difficult to make the initial distribution of GCs with a very small core. 
Moreover, simulations of a small r c also need higher spatial resolution than we currently 
employ, and so we leave this question for future study. 

4. Discussion and Conclusion 

Our main result is that steeper initial GC profiles experience stronger flattening as 
shown in Figure [5j a p ~ 3.5 is in boundary between strong and weak flattening. The 
results imply that if the initial slopes of both red and blue GCs are steeper than a p ~ 3.5, 
the difference in the slopes between two populations of GCs will become much smaller, 
independent of merger orbits. In particular, even only one dry merger can make the slope 
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flatter dramatically. Moreover, it also makes both slopes be around a r ~ 3.5 — 4.5. For 
example, if a Pire d = 5 and a Pi biue = 4, Figure [5] suggests that a rire d ~ 4.3 and a Pi bi ue ~ 3.6. 
The difference in the slope between red and blue GCs becomes significantly small. Therefore, 
the difference in the slopes of red and blue GCs can be a constraint of the number of major 
dry mergers. 

On the other hand, if the initial GC distribution has a slope shallower than a p ~ 3.5, 
the slope changes very little, and almost no change in the case of a p ~ 2. Therefore, if the 
initial distributions of red and blue GCs in merger progenitors are steeper than a p = 2, and 
the galaxies experienced a number of major dry mergers, it leads the distributions of both 
red and blue GCs to become close to a r ~ 2 for both red and blue GCs, and the difference 
in the slope of spatial distributions becomes difficult to be measured. 

It is also worth noting that if the initial core size is larger for blue GCs than for red 
GCs, but their initial slopes are the same, a major dry merger can make the slope for the 
blue GCs shallower than for the red GCs, as shown in Figure [H Therefore, the final slope is 
not a simple function of the initial slope. 



and blue GCs shallower and similar. Several studies ( 


s.e. iKissler-Patie 


1997; 


van den Bereh 


1998; 


Forbes 


2005; 


Lauer et al. 


2007; 


Emsellem et al. 


2007 


) claim that the various observed 



properties for ellipticals show a transition around M 



v 



-21, i.e. 



10 11 M . Some of these 



bimodalities have been interpreted in the picture of d ry merger hypo thesis for the growth o f 
massive ellipticals (e.g. ICapetti fc Balmaverdd 120061 ) . For example, lEmsellem et al.l (120071 ) 
show that slowly rotating ellipticals might be mainly affected by dry mergers. The slowly 
rotating ellipticals are more luminous than Mb ~ —20.5, while less luminous ellipticals are 
fast rotators. It may indicate that less luminous galaxies have not experienced any major 
dry mergers. Therefore, it is interesting to examine the observed distributions of red and 
blue GCs as a function of the stellar mass of ellipticals. If more luminous galaxies have 
experienced dry mergers, our results predict that the slopes for red and blue GCs in bright 
galaxies are shallower and similar. Unfortunately, the current observational samples are not 
enough to test it statistically. Below we provide some discussion based on the current limited 
measurements. 

NGC 1399 is one of e l liptic a ls whose spatial dist ributions of red and blue GCs are 
well-observed (IDirsch et al.l 120031 ; iBassino et al.l l2006al). When assuming M/L v = 5 and 
B — V = 0.9 for an elliptical galaxy, as used in iForbed (120051). the stellar mass of NGC 
1399 is about 5 x 10 11 M (M B = -21.8) farodie & Straderlhood ). The derived power-law 
indices of projecte d radial density profi les are 1.9 ± 0.06 and 1.6 ± 0.10 for red and blue 
GCs, respectively ( IBassino et al.ll2006al ). If the distributions are spherically symmetric, the 
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three-dimensional radial distributions have power-law slopes of a ~ 2.9 and « 2.6 for red 
and blue GCs, respectively. Note that the real slope might be slightly steeper, because 
the distribution of GCs is not infinite. The slope difference between red and blue GC 
distributions is Aa ~ 0.3 while both GC populations show a ~ 3. In addition to NGC 
1399, NGC 1407 is a brightest group gal axy with My = —21.86, and has a bimodal color 



distribution of GCs. iForbes et al.l (120061 ) report that the projected slopes of the red and 
blue GC density profiles are respectively 1.50 ± 0.05 and 1.65 ± 0.29, i.e., expected three 
dimensional slopes of a rec j = 2.50 at>i U e = 2.65, which are statistically the same as each other. 

On the other hand, NGC 1374 and NGC 1379 are less luminous than NGC 1399, having 
M y = —20.4 and —20.6 , respectively. The red and blue GCs of the two galaxies are studied 



in 



Bassino et al.l (l2006bl ). The derived projected density distributions of red and blue GCs in 
NGC 1374 have slopes of 3.2 and 2.3. Hence, the expected slopes in their three-dimensional 
density profiles are a re d ~ 4.2 and ctbiue ~ 3.3. The Aa in NGC 1379 is ~ 0.6, and the 
red and blue GCs have projected power-law slopes of ~ 2.9 and ~ 2.3, respectively. The 
low-luminosity galaxies have a larger difference in the density slopes of the red and blue GCs, 
and their slopes are steeper than the more luminous galaxies discussed above. Combining 
with our results, these four sample suggests that luminous galaxies, like NGC 1399 and NGC 
1407, may have experienced some dry mergers, while less luminous galaxies, such as NGC 
1374 and NGC 1379, may have not been formed through dry mergers, which is consistent 
with a scenario that the significance of dry mergers depends on the stellar mass. 

However, we als o find some giant galaxy that does not follow the above trend. Recently, 



Tamura et al.l (120061 ) measured the density profiles of red and blue GCs for giant ellipticals: 
M87 and NGC 4552. We have fitted the projected density profile in Figure 5 of their paper 
by a power-law profile, and the derived power-law index is 2.4 and 1.4 for red and blue GCs 
in M87 and 1.8 and 1.2 for red and blue GCs in NGC 4552, respectively. Although these 
giant ellipticals have relatively shallower slopes, the difference in slopes of red and blue GCs 
is significant. According to our result, these giant ellipticals are unlikely to have experienced 
significant dry mergers. Some giants might form without dry mergers, or some other factor, 
such as a number of minor mergers, might affect the slopes of GCs. 

Again, so far there are not many studies which derive the density distributions of both 
red and blue GCs. This problem limits the application of our results to the current obser- 
vational data. However, we expect that future observations will improve the statistics for 
nearby ellipticals, and our results will add an important framework to interpret the galaxies' 
evolutionary history. It will be also valuable to compare spatial distributions of blue and red 
GCs with other expected properties from dry mergers such as surface brightness profile. 
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Table 1. MERGER ORBITS 



Simulation 


Av (km/s) 


e 


r peri (kpc) 


A 


(206, 0.0, 0.0) 


1.0 


0.0 


B 


(205, 23, 0.0) 


1.0 


5.0 


C 


(193, 73, 0.0) 


1.0 


50 



Note. Eccentricity, e, and perigee dis- 

tance, r peri , are de rived from the reduced tw o- 
body problem (e.g. iKhochfar fc Burkertl 120061 ) 
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Fig. 1. — Initial number density (left) and velocity dispersion (right) profiles of GCs in 
merger progenitors. The input profiles of r c = 2.5 kpc and a p = 4.75 are shown as solid 
lines, while dashed lines represent the profiles of r c = 5.0 kpc and a p = 2.25. Dots are the 
profiles of sampled 10 4 GC particles that are well matched to analytical descriptions. 



-16- 



4 



3 



2 



1 1 1 1 1 1 1 1 1 


I | ▲ 


• a r c — c.OKpc 
oa r c =5.0kpc 


A *A O 


- 

- 








%/ 


- 










1 1 1 1 1 1 1 


1 I I I I 



2 3 4 5 

a 



Fig. 2. — Comparison of the initial slope («o) with the slope (ascyr) after 8 Gyrs of a 
simulation for an isolated galaxy without any interaction. Both high- resolution ( circle ) and 
low-resolution (triangle) cases show that the initial distributions of GCs are stable in a single 
isolated galaxy for 8 Gyrs. 
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Fig. 3. — Projected mass density distribution of stars in merger remnants. The de Vau- 
couleurs profile (line) is well matched to stellar distributions in all merger remnants. The 
distributions of the simulation A, B, and C are represented by triangle, rectangle, and circle, 
respectively. The measured effective radius is about 6.2 kpc. 
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Fig. 4. — GC distributions in merger remnants for simulation B. The initial distributions of 
r c = 2.5 kpc are shown for a p = 2.0,3.0,4.0, and 5.0 (solid line). The distributions of GCs 
in merger remnants are represented by dots. The initial distributions are scaled upward by 
a factor of two in order to match the total number of GCs in merger remnants. 
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Fig. 5. — a r (left) and a p — a r (right) for simulation B with r c = 2.5 kpc. The distributions 
in merger remnants are described as hqc oc r~ ar over 10 < r < 100 kpc. The dashed line 
represents no change in slope. 
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Fig. 6. — Slope of GC distributions in merger remnants produced from three different merger 
orbits. For the initial distribution at a fixed r c = 2.5 kpc, simulation C, with the highest 
angular momentum, shows the smallest change among the three orbits. 
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Fig. 7. — Distributions of GCs in three different merger remnants for a p = 5 and r c = 2.5 
kpc. The effect of different merger orbits is stronger in outer region than in inner region. The 
distribution in a merger progenitor is scaled upward by a factor of two for case of comparison 
(line). 



-22 - 




Fig. 8. — Comparison of the change of slopes (a p vs. a r ) between the high- resolution (square) 
and low-resolution run (circle) for simulation A with r c = 2.5 kpc. The change of slopes in 
a low-resolution run is comparable to those of a high-resolution run. 
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Fig. 9. — Dependence of GC distributions in merger remnants on the initial core radius, r c . 
Simulations B for r c = 2.5 and 5.0 kpc are compared. 
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Fig. 10. — Distributions of GCs in merger remnants for r c = 2.5 and 5.0 kpc. For simulation 
B and a p = 5, we present initial distributions in merger progenitors for r c = 2.5 (solid line) 
and 5.0 kpc (dashed line) and distributions in the merger remnant. 



